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Abstract 

We present a Monte Carlo study of the high-temperature phase of the two- 
dimensional driven lattice gas at infinite driving field. We define a finite-volume 
correlation length, verify that this definition has a good infinite-volume limit in- 
dependent of the lattice geometry, and study its finite-size-scaling behavior. The 
results for the correlation length are in good agreement with the predictions based 
on the field theory proposed by Janssen, Schmittmann, Leung, and Cardy. The 
theoretical predictions for the susceptibility and the magnetization are also well 
verified. We show that the transverse Binder parameter vanishes at the critical 
point in all dimensions d > 2 and discuss how such result should be expected in 
the theory of Janssen et al. in spite of the existence of a dangerously irrelevant 
operator. Our results confirm the Gaussian nature of the transverse excitations. 



1 Introduction 



At present, the statistical mechanics of systems in thermal equilibrium is quite well estab- 
lished. On the other hand, little is known in general for nonequilibrium systems, although 
some interesting results have been recently obtained jT]. It seems therefore worthwhile 
to study simple models which are out of thermal equilibrium. One of them was intro- 
duced at the beginning of the eighties by Katz, Lebowitz, and Spohn who studied 
the stationary state of a lattice gas under the action of an external drive. The model, 
hereafter called driven lattice gas (DLG), is a kinetic Ising model on a periodic domain 
with Kawasaki dynamics and biased jump rates. Although not in thermal equilibrium, 
the DLG has a time-independent stationary state and shows a finite-temperature phase 
transition, which is however different in nature from its equilibrium counterpart.^ 

Despite its simplicity, the DLG has not yet been solved exactly^ and at present there 
is still much debate on the nature of the phase transition pr7|ISlimiTmilHll2j . In Refs. jUlIZ] 
Janssen, Schmittmann, Leung, and Cardy (JSLC) developed a continuum theory which 
should capture the basic features of the transition and which provides exact predictions for 
the critical exponents. Several computer simulations in two and three dimensions [Tnill4[ 
IT^llTBj provided good support to these field-theoretical predictions, once it was understood 
that the highly anisotropic character of the transition required some kind of anisotropic 
finite-size scaling (FSS) [T7l[T3I. Still some discrepancies remained, prompting Garrido, 
de los Santos, and Munoz jHUHlCni to reanalyze the derivation of the field theory. On the 
basis of this analysis, they suggested that the DLG at infinite driving field should not 
behave as predicted by JSLC but should rather belong to the universality class of the 
randomly driven lattice gas (RDLG) fHllH]- This approach gives different predictions 
for the critical exponents that have been apparently verified numerically pm i21j. see also 
Ref. 122]. 

In view of these contradictory results, a new numerical investigation is necessary, 
in order to decide which theory describes the critical behavior of the DLG. For this 
purpose it is useful to consider quantities that are exactly predicted at least in one of 
the two theories. Here we shall focus on transverse fluctuations since in the JSLC theory 
transverse correlation functions are predicted to be Gaussian. Therefore, beside critical 
exponents, one can also exactly compute the FSS functions of several observables and 
make an unambiguous test of the JSLC theory. 

A basic ingredient of our FSS analysis is the finite-volume correlation length. In 
spite of the extensive numerical work, no direct studies of the correlation length have 
been done so far, essentially because it is not easy to define it. Indeed, in the high- 
temperature phase the model shows long-range correlations due to the violation of detailed 
balance [231211 • Therefore, no correlation length can be defined from the large-distance 
behavior of the two-point correlation function. A different definition is therefore necessary: 
In Refs. [221 cL parallel correlation length is defined. However, this definition suffers from 
many ambiguities (see the discussion in Ref. ^) and gives results for the exponent 

^For an extensive presentation of the DLG and of many generalizations, see Refs. [31^. 

^The DLG is soluble for infinite drive in the limit in which jumps in the direction of the field are 
infinitely more frequent than jumps in the orthogonal directions (H]- 
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which are not in agreement with the JSLC theory ^21- Even more difficuh appears the 
definition of a transverse correlation length because of the presence of negative correlations 
at large distances [T31 I3]. 

In this paper we define a finite-volume transverse correlation length generalizing the 
definition of the second-moment correlation length that is used in equilibrium systems. 
Because of the conserved dynamics, such a generalization requires some care. Here, we use 
the results of Ref. |26|. The main advantage of using a correlation length is the possibility 
of performing FSS checks without free parameters. We find that the susceptibility and 
the correlation length show a FSS behavior that is in full agreement with the idea that 
transverse fiuctuations are Gaussian. This immediately implies the JSLC predictions 
7j_ = 1 and z/_l = 1/2. We also study the FSS behavior of the magnetization, finding 
= 1.023(43), in agreement with mean-field behavior. Finally, we consider the 
transverse Binder parameter and find that it goes to zero at j3c as the volume increases. 
This supports again the Gaussian nature of the transverse mode and, as we discuss in 
detail, it is not in contradiction with the presence of a dangerously irrelevant operator. 
Indeed, since no zero mode is present in the theory, there should be no anomalous scaling. 

The paper is organized as follows. In Sec. El we describe the model and define the 
observables measured in the Monte Carlo simulation. Section El reviews the FSS theory 
and introduces the basic formulae that are used in the analysis of the numerical data. 
Next in Sec. El we consider the field theory of Refs. 011] in a finite geometry and compute 
the FSS functions for several observables. In Sec. |31 we describe the simulations and 
present the results which are then discussed and compared with recent findings pUl l^ 
in Sec. IHl We confirm the field-theoretical predictions of JSLC both for infinite-volume 
quantities and for the finite-size behavior. In App. ^we consider the 0{N) model for 
N oo above the upper critical dimension and discuss the critical behavior of different 
definitions of the Binder parameter. In particular, we show that, if it is defined in terms 
of correlation functions at nonzero momenta, then it vanishes for all d > 4: this is the 
same behavior as that expected in the DLG for ci > 2 on the basis of the JSLC theory and 
verified numerically for d = 2. In App. ^we sketch a one-loop calculation of the Binder 
parameter in the JSLC theory. A short account of the results presented in this work has 
been given in Ref. |27j . 

2 Definitions 
2.1 The model 

We consider a finite square lattice A and N particles, each of them occupying a different 
lattice site. A configuration of the system is specified by the set of occupation numbers 
of each site n = {rii G {0, l}}jgA- The standard lattice gas is characterized by a nearest- 
neighbor attractive ("ferromagnetic" in spin language) Hamiltonian 



where the sum runs over all lattice nearest neighbors. We consider a discrete-time 
Kawasaki dynamics j2HI, which preserves the total number of particles N or, equivalently, 




(1) 
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the density 

P^ = JX\^^'' (2) 

where |A| is the total number of sites in A. At each step, we randomly choose a lattice 
link If rii = rij, nothing happens. Otherwise, we propose a particle jump with 

probability w{AH/T), where 

AH = HA[n']-HA[n] (3) 

is the difference in energy between the new {n') and the old (n) configuration. If the 
probability w{x) satisfies 

w{-x) = w{x), (4) 

then the dynamics is reversible, i.e. satisfies detailed balance. Under these conditions 
there is a unique equilibrium measure given by 

Z-^{n'} ^ 

where (3 = 1/T. In the thermodynamic limit the lattice gas exhibits a second-order 
phase transition for pA = 1/2 and /3c = | In (v^ +1), which belongs to the standard Ising 
universality class. 

The DLG is a generalization of the lattice gas in which one introduces a uniform (in 
space and time) force field pointing along one of the axes of the lattice, i.e. E = E^t: It 
favors (respectively suppresses) the jumps of the particles in the positive (resp. negative) 
x-direction. If A is bounded by rigid walls, then E is a conservative field and it can be 
accounted for by adding a potential term to HA[n\. Therefore, the system remains in 
thermal equilibrium. The net effect of E is simply to induce a concentration gradient in 
the equilibrium state. 

Here, we consider instead periodic boundary conditions.^ In this case, the field E does 
not have a global potential and the system reaches a stationary state which, however, is 
not a state in thermal equilibrium. 

In the DLG transition probabilities take into account the work done by the field during 
the particle jump from one site to one of its nearest neighbors. In this case one proposes 
a particle jump with probability w{[3 AH + PE£) , with i = (—1, 0, 1) for jumps (along, 
transverse, opposite) to ic. 

For E and p\ = 1/2, the system undergoes a continuous phase transition at 
an inverse temperature Pc{E) which saturates, for E oo, at /?c(oo) ^ 0.71/5c(0). For 
P < Pc{E) particles are homogeneously distributed in space, while for /? > Pc{E) phase 
separation occurs: Two regions are formed, one almost full and the other one almost 
empty, with interfaces parallel to E. 

■^In principle, it is enough to consider periodic boundary conditions in the field direction. The boundary 
conditions in the transverse directions are largely irrelevant for the problems discussed here. 
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2.2 Observables 

We consider a finite square lattice of size L\\ x L\_ witli periodic boundary conditions. We 
define a "spin" variable sj = 2nj — 1 and its Fourier transform 

0(k)^5^e^'^-J.j, (6) 

jeA 



where the allowed momenta are 



_ / 27rn 27rm\ 

with {n,m) G Z^^n x Z^.^^. 

We consider the model at half filling, i.e. for = 1/2. Then 

5^Sj=0, i.e. 0(ko,o)=O. (8) 

In the ordered phase |0(k)| takes its maximum for k = ko.i, and the expectation value on 
the steady state of its module 

m(/?;L||,L^)^i^(|0(ko,i)|) (9) 

is a good order parameter. 

In momentum space the static structure factor, the Fourier transform of the two-point 
correlation function, 

G'(k;L||,L^) = i^(|0(k)r) (10) 

vanishes at ko,o because of Eq. (jHJ and attains its maximum at kg^i, so that it is natural 
to define the susceptibility as'' 

X±(/?;L||,L^)=G(ko,i;L||,L^). (11) 

We also define the four-point connected correlation function 

G'W(ki,k2,k3,k4;L||,L^) = i^(0(ki)0(k2)0(k3)0(k4))co„n , (12) 

and the related transverse Binder cumulant g{(3; L\\, L±) defined as 

„r« r r mKy)\'} _ GW(k.,.. ko,i. -kp,,. -k„,i; L,, L^) 

''"'•^"■^-) = 2-(Wk„,0|T" |A||G(k„,,iL,|.L,p ■ 



^We note that the susceptibihty defined by using the hnear response theory does not coincide in 
nonequihbrium systems with that defined in terms of the Fourier transform of the two-point correlation 
function. 
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Next, we would like to define a correlation length. In infinite- volume equilibrium systems 
there are essentially two different ways of doing it. One can define the correlation length 
in terms of the large-distance behavior of the two-point correlation function or by using 
its small-momentum behavior (second- moment correlation length). In the DLG the first 
method does not work. Indeed, in the high-temperature phase the two-point correlation 
function always decays algebraically with the distance. Moreover, it is not positive definite 
because of negative correlations in the transverse directions j3j . This peculiar behavior is 
due to the fact that in the infinite- volume limit (at fixed temperature) the static structure 
factor G(k; oo, oo) has a finite discontinuity at k = 0. 

In this paper we propose a new definition that generalizes the second-moment corre- 
lation length used in equilibrium spin systems. The basic observation is that in the DLG 
the infinite-volume wall-wall correlation function decays exponentially, i.e. 

^G((a;||,xx);oo,oo) = /" d'^-^g^ (^((0, g^); oo, oo)e*'?^-^^ ~ e""!^^!, (14) 

as in equilibrium systems. This holds at tree level both in the JSLC and in RDLG field 
theories and to all orders of perturbation theory in the JSLC theory, see Sec. El Therefore, 
we can define a correlation length as in equilibrium systems, paying due attention to the 
conserved dynamics. Here, we follow Ref. [20] , where we discussed the possible definitions 
of correlation length in the absence of the zero mode, as it is the case here. 

We consider the structure factor in finite volume at zero longitudinal momenta 

G^{q-L\\,L^)= G((0,g);L||,L^), (15) 

(note that the conservation law implies G±{0] L\\, L±) = 0) and define a finite- volume 
(transverse) correlation length^ 



Gxte;L||,Li) J 

where g„ = 2sin(7rn/Lj_) is the lattice momentum. Since G±{0; L\\, L±) = 0, qi and qj 
must not vanish. Moreover, as discussed in Ref. |^, the definition should be valid for 
all (3 in finite volume. Since the system orders in an even number of stripes, for i even 
G±{<li', L\\, L±) is zero as /3 ^ oo. Therefore, if our definition should capture the nature 
of the phase transition, we must require i and j to be odd. Although any choice of i,j is 
conceptually good, finite-size corrections increase with a phenomenon which should 
be expected since the critical modes correspond to q —>■ 0. Thus, we choose = (1, 3), 
defining = ^13. 

^In Ref. |l26j we showed that any good finite- volume correlation length must satisfy two properties: (i) 
it must be finite for all T 7^ and L±, Ly < 00; (ii) it must diverge as T — > even in finite volume. The 
definition (|16|l satisfies these two properties. Note that in infinite volume one can also define a correlation 
length from the large-distance behavior of the correlation function, i.e. one can define S,od.± — where 
K is defined in Eq. p4|l . Such a definition is not convenient here, since it does not admit a finite- volume 
generalization. 



^ij{L\\,L_L) 



\ 



Qi - Qi 
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The same method can be used to define a longitudinal correlation length, although in 
this case we do not have an all-order proof (not even in the JSLC theory) that wall-wall 
longitudinal correlations decay exponentially. It is enough to consider the longitudinal 
structure factor at zero transverse momentum 

G||(g;L||,Lx) = G((g,0);L||,Lx), (17) 

and use again Eq. p6|) . In this case, there is no reason to avoid the use of even i or j and 
thus we define = ^12- 

Near a phase-transition point the quantities we have defined above show power-law 
divergences. As usual, see, e.g., Ref. we define transverse exponents by assuming 

X± ~ t-^\ (18) 
for t = {(3c — l3)/l3c —>■ 0^. The magnetization vanishes in the low-temperature phase as 

mr^{-tf^. (19) 



3 Finite-size scaling 

In the neighborhood of a critical point the behavior of long-range observables is controlled 
by few quantities, corresponding in the renormalization-group language to coordinates 
parametrizing the relevant directions in the infinite-dimensional coupling space. When 
the system is finite, its size plays the role of another relevant operator. This means that 
an observable O, which diverges in the thermodynamic limit as 

Ooo(/9)~t-^° fort=l-|-^0+, (20) 

behaves in a finite system of size Ly x 

0{(3; L\\,L^) ^ r^^hoit'VLr, S) ^ L^f''f2,o{t-''/L^-, S) ^ L^f'''h,o{U{P)/Lu S), 

(21) 

where ^00 (/5) is the infinite- volume correlation length and S = L\\/L±_ is the aspect ratio 
that is kept fixed in the FSS limit. From Eq. 1)211) we can derive a general relation for the 
ratio of 0{f3; L\\, L±) at two different sizes (L||,L_l) and {aL\\,aL±). In the FSS limit we 
obtain 



0(13; L,,L^) 

where we have replaced ^oo{P)/Lj_ with ^{P; L\\, L±)/ L± by inverting ^{P;L\\,L±) ^ 
L±f3,^{^oo{P)/Lj_; S). The function Fo{a,z,S) is universal and is directly accessible nu- 
merically, e.g., by Monte Carlo simulations — no need to fix any parameter — since all quan- 
tities appearing in Eq. ()22|) are directly measurable. Moreover, as we shall show below, all 
critical exponents can be determined from the FSS functions FQ{a,z,S) independently 
of the critical temperature. 
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If we define z = ^{P; L\\, L±)/ L±, then z varies between and z*{S), where z*{S) is 
defined by 

z*iS) = fs,doo,S), (23) 

or imphcitly from 

a = F^{a,z*{S),S). (24) 

The value z*{S) is directly related to the behavior of the finite-size correlation length at 
the critical point, since ^{Pc, L\\, L±) ^ z*{S)L±. For finite-temperature phase transitions 
z*{S) is finite. By considering the behavior of the FSS functions at z*{S) we can determine 
the exponents 70 / 1^- Indeed, at the critical point we have 



0(/3e;L||,L^)~LI°/^ (25) 



so that 



and therefore 



Foia, z iS), S) = = , (26) 



70 logFo(«,^*(^),^) 



(27) 



u log a 

By studying the behavior of F^{a, z, S) in a neighborhood of z*{S) it is also possible to 
derive the exponent u. Using the fact that 



near the critical point, we obtain 

dF^{a, z, S) 



z*{S) + a{S){P-P,)LY'', {21 



dz 



= a{a^/'' - 1). (29) 

Z=Z*{S) 



The above-presented results are valid for an isotropic system. On the other hand, the 
numerical simulations and the field-theoretical studies predict that the phase transition 
in the DLG is strongly anisotropic. For example, the scaling form of the critical static 
two-point function should be 

G{k\\, k^) ^ ii'^+''^G{ii^+%\, fik^), (30) 

where rj± is the anomalous dimension of the density field (see Ref. |l3| for definitions) and 
A is the so-called anisotropy exponent. 

It is then natural to assume the existence of two correlation lengths which diverge 
with different exponents i'± and i/y related by |SI 

= (1 + A)z/x. (31) 

These considerations call for an extension of the FSS arguments. A phenomenological 
approach to FSS for the DLG has been developed [221, keeping into account the strong 
anisotropy observed in the transition (for d = 2 see Refs. [T3|ll4pi7j . for d = 3 see Ref. ^16j). 
Following this approach, we assume that all observables have a finite FSS limit for Ly, 

— i> 00 keeping constant: 
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• the anisotropic aspect ratio S'a = L^^^^^~^^^ /L±; 

• the FSS parameter ^±,oo{P)/ L± (or equivalently its longitudinal counterpart). 

Then, Eq. pT|l still holds by using the correct parameters, i.e. by replacing S with Sa, 
V with and ^oo with ^_l,oo- Analogously Eq. (j22I) is recast in the form 

0(/3;L||,Lx) 

Equation p2|l is the basis for our analysis of the phase transition in the DLG . For the 
transverse finite-volume correlation length we use ^13 defined in Sec. 12.21 

In anisotropic systems the FSS limit must be taken at fixed S'a, which in turn requires 
the knowledge of the exact value of A. It is thus important to understand how we can 
single out the correct value of A from the simulations. First, it should be noticed that 
observing FSS does not imply that we are using the correct value of A jSIj. Indeed, note 
first that Eqs. (|2ip and (j32|) still hold for L\\ = 00 — this corresponds to S'a = 00 — i.e. for 
a geometry 00 x (in two dimensions this is a strip). Then, imagine that we use an 
incorrect value 6, keeping 5*^ fixed in the FSS limit. If 5 > A, Ly increases much faster 
than it should and 00. We thus expect to obtain the FSS behavior corresponding 

to a geometry 00 x L^~^. Thus, we should be able to observe scaling whenever we use L± 
as reference length and consider transverse quantities. On the other hand, longitudinal 
quantities are expected not to scale properly. For instance, should not have a good 

FSS behavior when plotted vs ^±/L±. Indeed, at fixed Sa we have that ~ Ly ~ L^^^. 
If instead Ss is fixed and 6 > A, L\\ increases too fast, so that we expect to be controlled 
by the transverse size only. Thus, ^y should increase slower than Ly, and thus, at fixed 
one should observe ^||/-^|| — >0. If5<A the same argument holds by simply 
interchanging longitudinal and transverse quantities. This observation provides therefore 
a method to determine the correct value of A. It is the value for which both correlation 
lengths scale correctly, i.e. ^\\ ~ Ly and C,± ^ L±. 




4 Field-theory description of the DLG 

As we explained in Sec. 12.11 the DLG is a lattice gas model. However, in a neighborhood 
of the critical point (critical region) we can limit ourselves to consider slowly-varying (in 
space and time) observables. At criticality the lattice spacing is negligible compared to 
the length and time scales at which long-range order is established so that it is possible to 
formulate a description of the system in terms of mesoscopic variables. In principle, the 
dynamics of such variables can be obtained by coarse graining the microscopic system. 
However, given the difficulty of performing a rigorous coarse-graining procedure, one 
postulates a continuum field theory that possesses all the symmetries of the microscopic 
lattice model. By universality the continuum theory should have the same critical behavior 
of the microscopic one. 

Unfortunately, there is at present no consensus on the field theory that describes the 
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critical behavior of the DLG.^ The theory originally proposed by JSLC [7lin] has been 
recently disputed in Refs. [21131101 (see also Ref. [H]), where it is proposed that the DLG 
at infinite driving field is in the same universality class of the randomly driven lattice gas 
(RDLG) ^HlEl- Both theories agree on the anisotropic nature of the phase transition, 
but disagree on its origin and make different predictions for the universal quantities. 

In the JSLC theory the critical exponents are exactly predicted, and, for 2 < d < 5, 
they are given by 

V± = 0, (33) 

= (34) 

7± = 1, (35) 

P± = \. (36) 

A = \{S-d). (37) 

In the RDLG model critical exponents are only known up to two loops perturbatively in 
e = ?) — d. Therefore, it is difficult to estimate how much they differ from the JSLC ones. 
We note that A should be quite different in the two theories in two dimensions. In the 
RDLG we should have A ^ 1, since jT2] A = 1 — 77/2 and we expect 77 to be small. On 
the other hand, Eq. gives A = 2. 

For the JSLC theory, not only do we have exact predictions for the exponents, but 
we can also compute exactly the transverse structure factor (ITHjl . Keeping into account 
causality [HIESIEIIES] and the form of the interaction vertex one can see that for k\\ = 
there are no loop contributions to the two-point function G{k) (and also to the two-point 
response function). Thus, for all 2 < < 5, G±{k) is simply given, in the field-theoretical 
approach of JSLC, by the tree-level expression 

G^ik) = (38) 

where r is a squared "bare mass" that vanishes at criticality. Two observations are in 
order. First, r = bt + Oit^) for t = (/^^ — /?)//?c — ^ with h positive constant. Second, 
the function that appears in Eq. (j!?H|) refers to the coarse-grained fields, which, in the 
critical limit, differ by a finite renormalization from the lattice ones. Thus, for the lattice 
function we are interested in, in the scaling limit t — 0, — > 0, with k'^ /t fixed, we have 

Gxa«(*) = (39) 

where Z and h are positive constants. 

On the same footing, we can conclude that all correlation functions with vanishing 
longitudinal momenta behave as in a free theory. In particular, the Binder cumulant 



'The effects of an external drive on the standard Model-B dynamics were also studied in Ref. |32| . 
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defined in Eq. fll3|) vanishes. It is also important to notice that Eq. (jH^ imphes the 
exponential decay of 



G^M^±) = j d^-'ke'^-^ G^,iatt(A:), (40) 

which fully justifies our definition of transverse correlation length. 

Field-theoretical methods can also be used to determine the FSS behavior of the 
model. Predictions are easily obtained by following the method applied in equilibrium 
spin systems (see, e.g., Ref. jHSl Chap. 36] and references therein). The idea is quite 
simple. Consider the system in a finite box with periodic boundary conditions. The finite 
geometry has the only effect of quantizing the momenta. Thus, the perturbative finite- 
volume correlation functions are obtained by replacing momentum integrals by lattice 
sums. Ultraviolet divergences are not affected by the presence of the box and thus 
one can use the infinite-volume renormalization constants. Once the renormalization is 
carried out, one obtains the geometry-dependent finite-size correlation functions. 

The considerations we have presented above for the infinite- volume case apply also in 
finite volume and thus Eq. (|39p holds in this case. Using Eq. (|39|) . in the FSS limit we 
find 

^^^^V^"'^^^ = [{27r? + btLl]-'/\ (41) 

valid for t — 0, Ly, L±_ — > oo with tL^ fixed. Although we have not explicitly mentioned 
Sa and such a quantity does not appear in Eq. (jH)), this expression is expected to be 
valid only if the FSS limit is taken by keeping S'a constant. Taking the limit L± oo 
we obtain ^±,oo(/^)^ ~ ^/{bt)> so that we can write Eq. (|41|) also in the form 

1 _ 1 47r2 

[UP;L\\,L^W ~ ^ ^ ■ ^^^^ 

Using Eqs. ()42|1 and (p^ we can also compute the scaling functions F^{a, z, Sa) and 
F^{a, z, S/^) defined in Eq. (jH^ . We obtain 

F^{a,z,SA) = [l-(l-a-')(27r)V]"'/', (43) 
F^(a, z, Sa) = F^{a, z, Sa) = [l - (l - a'^) {2tt)^z^] . (44) 
Note that these functions do not depend on Sa and that 

z*{Sa) = ^ . (45) 

A peculiarity of the JSLC theory is the presence of an operator with renormalization- 
group dimension 2a = 2{d — 2)/3 that is dangerously irrelevant for 2 < < 5 and 
becomes marginal at = 2. Keeping into account the coupling u associated with the 
dangerously irrelevant operator, in the JSLC theory we have [Hj 



rs,n({(gii,g±)},c^,i-,M;^ii,^±) = (46) 
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where is the one-particle irreducible correlation function of n density fields and n 
response fields. For the transverse structure factor it implies 

G^{q^, T, u- , L^, = fG^iiq^, fr, ^^'^u- r'L^, r^-^Lu). (47) 

Setting i = L± and expanding for L_l oo we obtain for 2 < < 5 

r, u; , L^, Ly) = Llh{q^L^, tLI; Sa)[1 + ^(nL-^-)]. (48) 

Now the leading term is given by Eq. (j38p . and thus, in the absence of zero mode, which im- 
plies |g±|Iv_L > 2n, f2{q±Lj_,TL\; Sa) is regular and finite in the whole high-temperature 
phase. Thus, the dangerously irrelevant coupling can be neglected for all 2 < < 5. 
Note that this argument would not apply to the zero mode if it were present. Indeed, for 
q± = 0, f2{0,TL\] Sa) = '^/tL\ which is singular as r — *• 0, giving rise to an anomalous 
behavior. In App. ^we discuss this phenomenon in the large- limit of the 0{N) model 
above the upper critical dimension, i.e. for d > 4, showing that no anomalous behavior 
is observed in the absence of zero mode (in particular the Binder parameter vanishes).^ 
Therefore, for d > 2, Eqs. PHj) and ()44|) should hold without changes. For d = 2 the 
operator becomes marginal {a = 0) and therefore we expect logarithmic corrections to 
the formulae previously computed. In the absence of any prediction, we will neglect these 
logarithmic violations. As it has been observed in previous numerical studies, if present, 
they are small ^3]. As we will discuss, this is confirmed by our numerical results. 

Finally, we wish to compute the behavior of the Binder parameter keeping into account 
the presence of the dangerously irrelevant operator. Considering the connected static four- 
point correlation function for vanishing parallel momenta, Eq. (j46|) implies 

G(^)({(0, gx)}, r, u; L\^,L^) = i''+^+^G^'\{{0, iq^)}, i\ r^'^u; r'-^L\\,r'L^). (49) 

Now is of order u, and therefore, setting i = L±, q± = ±q with |g| = 2n/L± as in 
the definition of the Binder parameter, cf. Eq. (|13|). we obtain 

G(^)({(0, ±g)}, r, u; L\\,L^) = uL\U{tLI, Sa)[1 + 0{uLl^-)]. (50) 

In App. ^we have verified this equation to first order in perturbation theory, proving 
also that /4(0, Sa) 7^ 0. It follows that the Binder parameter behaves as 

r5-d 

^7(/5;^ll,^±) = -^h{TLl,SA){^7,^ + TLlf[i + 0{uLi'")] 

= uLl^'^f,{TLl,SA)[l + 0{uLl^'^)l (51) 

with /g(0, Sa) 7^ 0. Therefore, for all 2 < d < 5 the Binder parameter vanishes, in spite 
of the presence of the dangerously irrelevant operator. In (i = 2 we expect logarithmic 
corrections. Note again that g{[3c] L\\^, L^) vanishes as L_l — >• cx) because /2(0,S'a) 7^ 0, a 

^Even in the absence of zero mode it is still possible to observe anomalous scaling by performing a 
noncanonical scaling limit. One should consider r — > 0, L± with tL\ + Aii^ ^ at the same time. 
See Ref. for a discussion in the A^-vector model. 
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consequence of the fact that the Binder parameter we use here is defined at nonvanishing 
momenta. 

Finally, we wish to discuss the distribution function of the order parameter = 0(ko_i) 
in the JSLC theory. For [3 ^ [3c and Ly, L_l oo, such a quantity has a Gaussian 
distribution, i.e. 

iVexp [- 1^1' ) dijdr. (52) 

where is a normalization factor, and indeed the Binder parameter vanishes in this 
limit. Since the dangerously irrelevant operator is truly irrelevant, such a distribution 
is also valid at the critical point and in the FSS limit, allowing the computation of the 
distribution function of the magnetization. We obtain 

= -r 3— r, (5o) 

4L||L;[-^ ^ ' 

with corrections of order L'^'^ (logarithms in c? = 2). Then, we obtain 

F^(a, z, Sa) = ^(rf+A)/2 -^x(^> SaY^^, (54) 
with logarithmic corrections in d = 2. 



5 Numerical Simulation 
5.1 Setup 

In order to study the critical behavior of the DLG in two dimensions, we perform an exten- 
sive Monte Carlo simulation. We use the dynamics described in Sec. 12. II with Metropolis 
rates, i.e. we set 

w{x) = min (1,6""^). (55) 

Simulations are performed at infinite driving field. Therefore, forward (backward) jumps 
in the direction of the field are always accepted (rejected). 

The dynamics of the DLG is diffusive and the dynamic critical exponent associated 
with transverse fluctuations, which represent the slowest modes of the system, is expected 
to be close to 4 both in the JSLC and in the RDLG models. In the JSLC theory jS] 
z± = 4 exactly . Thus, it is important to have an efficient implementation of the Monte 
Carlo algorithm in order to cope with the severe critical slowing down. 

We use a multi-spin coding technique, evolving simultaneously A^muiti independent 
configurations. We took particular care in optimizing the value iVmuiti- On Pentium and 
PowerPC processors — the computers we used in our simulations — we observed a non- 
monotonic behavior of the speed with A^muiti- For instance, on a Pentium processor, we 
reach a speed of 1.3 ■ 10® spin-fiips/sec for A^muiti = 32 and 2.7 ■ 10® for iVmuiti = 128. The 
speed drops to 2.0 ■ 10® for A^muiti = 192 and then increases again, reaching 3.8 ■ 10® for 
^muiti = 960. However, by increasing A^muiti we increase the memory and disk requirements 
so that we used A^muiti = 128 as a good compromise. 
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For the pseudo-random numbers we use the Parisi-Rapuano congruential generator, a 
32 bit-shift-register generator based on 

O-n = (O'n-24 + ^n-Ss) X0Ra„_61, (56) 

where a„ is an unsigned 32-bit integer. 

The main purpose of our work is to test the theoretical predictions of the JSLC field- 
theoretical model discussed in Sec. 01 It predicts A = 2 and therefore we have performed 
simulations on a set of lattices with L\\/L\ constant. We considered the following val- 
ues of (Ly, L^): (21,14), (32,16), (46,18), (64,20), (88,22), (110,24), (168, 28), (216, 30), 
(262,32), (373,36), (512,40), (592,42), (681,44), (778,46), (884,48). It is easy to verify 
that ^2 ^ 0.200 (more precisely 0.197 < ^2 < 0.202 for < 28 and 0.19992 < ^2 < 1/5 
for > 30). We considered several values of p: 0.28, 0.29, 0.3, 0.305, 0.3075, 0.31, 
0.3105, 0.311, 0.31125, 0.3115, 0.31175 0.312. As we shall see below, all lie in the disor- 
dered phase, albeit very near to the critical point. For a few values of (3 we have also 
performed simulations for different values of {L\\,L±). We have considered a sequence of 
lattices with ^2 ~ 0.300 and ^2 ~ 0.100 [the largest lattice has (L||,Lx) = (1350,30), 
(414,48) respectively] and a sequence with A = 1 and 5*1 ~ 0.106 [the largest lattice has 
(L||, L±) = (245,48)]. The results for A = 2 and 5*2 ~ 0.200 are presented in Tables [T] and 

121 

It is very important to be sure that the system has reached the steady-state distribution 
before sampling. Metastable configurations in which the system is trapped for times much 
longer than typical relaxation times in the steady state are a dangerous source of bias. 
In the DLG, configurations with multiple stripes aligned with the external field are very 
long-lived and may persist for times of the order of a typical simulation run. To avoid 
them, we started the simulations for the largest systems from suitably rescaled thermalized 
configurations of smaller lattices at the same temperature and value of S'a- 

We have performed a detailed study of the dynamic correlations. For P = 0.312, 
the value of (3 that is nearest to (3c, we compute the autocorrelation time for the 
susceptibility x±, which is expected to have a significant overlap with the slowest modes 
of the system. The results are reported in TableEl The autocorrelation times are expressed 
in sweeps, where a sweep is conventionally defined as the number of proposed moves equal 
to the volume of the lattice. It is easy to verify that ~ L\_ as expected. For each (3 
and lattice size we always make at least 3-10'' sweeps, so that our runs are at least, but 
in most of the cases much longer than, approximately 10^ r^^-. 

The statistical variance of the observables is estimated by using the jackknife method 
To take into account the correlations of the samples, we used a blocking technique in the 
jackknife analysis, using blocks of 256 x 10^ sweeps. They are much longer than typical 
autocorrelation times and thus different blocks are statistically independent. 

5.2 Results 

First, we check that the correlation length ^_l(/5; Ly, L^) has a good thermodynamic limit, 
independent of the aspect ratio, i.e. that, for j3 < jSc, ^±(/5; L\\, L^) has a finite limit when 
L||, L± — > oo in an arbitrary way. 
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Instead of ^±{f3; L\\, L_i_) we consider (the reason of this definition will become clear 
below) 

Of course, if ^±{P; L\\, L±) has a good thermodynamic limit, also r(/3; Ly, L±) has a finite 
limit and 

hm r(/3;L||,Lx) = -^^. (58) 

In Fig-^we show r(/3; Ly, L_|_) for four values of (3 and two different sequences of lattice 
sizes, one with 5*2 ~ 0.200 and the other one with 5*1 ~ 0.106. The numerical results show 
that t{(3; L\\, L^) has a finite infinite- volume limit, which is reached from below using 
the data with fixed 5*2 and from above using the data with fixed Si. The corrections 
to Eq. are not clear. In equilibrium systems the second-moment correlation length 
converges to ^±,oo(/5) with corrections that decrease as 1/-Z^y, 1/-^!) although in many 
cases, with an appropriate finite-volume definition, such corrections are so small that an 
apparent exponential convergence is observed jSHl- These corrections can be easily related 
to the behavior of G'_l,oo(?) in infinite volume. Indeed, Gj_{q] L\\, L±) converges to G±,oo(q') 
exponentially and thus we obtain for ^± = ^13 defined in Eq. flTB|) 



ei(/3;L||,L^) = eloo 



^loo / ,_2 , IOtt^ ^ 2071^C{(3) 



^1 + ^ |^-4vr^ + + ]+0{L- 

where 

c{f3) 



(59) 



(60) 



q=Q 

This expression is valid for any f3, not only near the critical point. It shows that corrections 
vanish as Lj^, although it does not allow to compute them, since c(/3), which for scaling 
reasons is expected to scale as as /3 — > Pc, is unknown. In the DLG case, we do 
not know what are the finite-volume corrections to the structure factor, but it is still 
reasonable to conjecture that corrections to C,± vanish as L^^. The results reported in 
Fig- E confirm this behavior. One apparently observes a Lj^ correction, especially for the 
lattices with fixed 5*1. It is interesting to observe that, considering the data with fixed 
5*2, the finite-size corrections decrease as /5 — > jSc- This is due to our specific definition of 
t{(3; L\\, L±) and is the motivation for this choice. Indeed, Eq. (jl^ implies 

r{P;L\\,L^)^^^^+oiL~/), (61) 

i.e. corrections vanish faster than Lj^, in the FSS limit at fixed 5*2. The data in Fig. H 
confirm this behavior and thus support the JSLC prediction (jl2I). 

For (3 = 0.311 we have performed a more detailed check by comparing results for three 
different values of 5*2. We define 

6{P- L\\,L^) = r(/5; L,,, L^) - t^{P), (62) 
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where Toq{P) is the extrapolated value of t{P; L\\, L±) for L±_ = oo. If Eq. holds, 
such a quantity vanishes as Lj^"*^ for L± ^ oo, where a; > is a correction-to-scaling 
exponent. In Fig.|21we report LI^Sl^P; L\\, L^) for three different values of a. For a = 2+uj, 
points with constant 5*2 should lie on an approximately horizontal line. From Fig. |21 we 
obtain u ^ 1, i.e. r(/5; L\\,L±) = Too{i3) + 0(L^^) for (3 close to jSc- 

We wish now to perform a detailed FSS analysis comparing the numerical data with 
the predictions PHj) and ()44|) . First, we consider the correlation length using a = 2. 
In Fig. ISlwe report our results for ^±{P;8L\\,2L±)/^±{P] L\\, L±) vs. ^±{P] L\\, L±)/ L±. 
The solid line is the theoretical prediction (j43|) . It is clear that, as the size of the lattice 
increases, the points converge towards the theoretical line. We would like to emphasize 
that in this plot there are no tunable parameters. Thus, the observed collapse is very 
remarkable. To get rid of the small corrections to FSS that are still present, in Fig. |3] 
we present the same data using on the horizontal axis ^±{P;SL\\,2L±)/ {2L±). By using 
the values of ^± corresponding to the larger lattice, scaling corrections are systematically 
reduced. Note that the very good agreement between theory and numerical data gives 
h'± = 1/2. Indeed, using Eq. (jIHj) and Eq. we obtain the mean- field value for the 
exponent z/^. 

Next, we checked the FSS behavior of the susceptibility. In Fig. we report our 
numerical data together with the theoretical prediction (j44j) . Again, we observe a very 
good agreement. This result immediately implies 7_l = 1, cf. Eq. (j^7j) . 

In Fig. ini we present the same plot for the magnetization. The results are in rea- 
sonable agreement with the theoretical prediction ()54|). although in the region of small 
^(/?; L_l)/L_l one can see some deviations. As we discuss below, such deviations can 
be interpreted as corrections to scaling, that in our case decay quite slowly, probably 
as an inverse power of logL^. If we perform a polynomial fit of the data, we obtain 
Fm{2, z\ S2) = 0.491(15) for z* = 1/{2tt). Using Eq. ^ we obtain 

— = 1.023(43). (63) 

Such a result is in very good agreement with the JSLC prediction P^/u^ = 1. 

The results for the Binder parameter reported in Fig. [7| show a good scaling behavior, 
falling onto a single curve with small corrections. By using Eq. we obtain approx- 
imately g{f3c', L\\, L±) ~ L^''^^^^^\ Thus, g = a.t the critical point, in agreement with 
the JSLC theory that predicts Gaussian transverse fluctuations at the critical point. The 
value of the exponent is difficult to interpret and we suspect that the Binder parameter 
is going to zero as some negative power of logL^, which however cannot be distinguished 
numerically from a small exponent. 

Finally, we report results for the parallel correlation length. Fig. |H1 reports a plot of 
(^ll(/?; L||, L^)/L|| vs. ^j_(/5; L||, L^)/L_|_. There are very large corrections to FSS but the 
data apparently collapse on a well-defined curve as the size of the lattice increases. This 
is again a very important test of the JSLC theory. Indeed, as discussed in Sec. IHlthe data 
for the correlation length can both scale only if we have correctly chosen the anisotropy 
exponent A. Thus, Fig. |Hl obtained using data with A = 2, strongly supports such a 
value for the anisotropy exponent. 
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In the JSLC theory the order parameter has the Gaussian distribution ()52|). However, 
such an expression does not take into account the presence of the leading corrections and 
therefore the FSS behavior of the Binder parameter. If we include these terms we expect 
a more complex form 

Nexp[-V{\tlj\'^)]diljdij*. (64) 

For small it is sensible to expand V'd'i/'p) in powers of lipl"^ and thus, if only the 
lowest moments of {ipl"^ are of interest, we can try to approximate Vd'i/'p) with its first 
two terms, 

V{\ij\^)^a\^\^ + b\^\\ (65) 

where b/a^ ~ uL'^'^ times a function of tL\. 

We now show numerically that this approximation describes remarkably well our re- 
sults for m, xx, and g: Essentially all corrections to scaling we observe can be taken into 
account by simply assuming the distribution function and (jHH)- 

If we define 

Mn{z) = ■' , , ' 1,12 1,14 , 66 

then the approximation predicts 



where z = h/a?. Thus, X turns out to be a function of the Binder parameter g, implicitly 
defined by Eq. (jFTfj) . In Fig.lHlwe report X versus g together with the theoretical prediction 
(I67|) . The agreement is remarkable, indicating that all corrections are taken into account 
by a simple generalization of the distribution of the order parameter. Note that, as (7 — *• 0, 
X approaches 7r/4 ^ 0.785, which is the value expected for a purely Gaussian distribution. 

Finally, we compute (3c from the critical behavior of ^±,oo(/3). In order to compute the 
infinite- volume correlation length, we can use two strategies. One consists in extrapolating 
r(/3;L||,L^) to L± — > cxd at fixed 82- To minimize corrections to scaling we only use 
the data with ^2 ~ 0.200. In Table HI we show the results of the fits r(/3; L|[, Lj_) = 
Too(/?) + <i{l3)L^ for each (3. In all cases we performed several fits, including each time 
only the results with L > Lmin for increasing Lmin. The results we report correspond to 
the fit for which the is reasonable, i.e. it is less than the value corresponding to a 
95% confidence level. In two cases we have not been able to fulfil this criterion, with a 

which is however only slightly larger. The value of Tao{(3) gives the estimate of the 



infinite-volume correlation length ^±,oo(/5) = 1/vToo(/3). 

In the previous analysis we have not made any assumption. More precise estimates of 
^±,oo(/5) can be obtained if we assume the validity of the JSLC theory and use the method 
of Ref. to determine infinite- volume quantities.® Such a method is particularly efficient 

similar method has been introduced in Refs. |41ll42j . 
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and it has been successfully applied to many different equilibrium models 
ElEHlllS]- For the calculation we assume here the theoretical prediction for the function 
F^{a, z, S^), cf. Eq fliH|) . and use the iterative algorithm with a = 2. To avoid scaling 
corrections we only consider the data with L > L^^[^ = 30 and 5*2 ~ 0.200. Then, for 
each P and L± we obtain estimates of C,_l^oo{P)- Results corresponding to the same /3 are 
then combined to give the final estimate reported in Table El Here N is the number of 
degrees of freedom (the number of data with the same (3 minus one) and is the sum of 
square residuals which gives an indication of the consistency of the various extrapolations. 
The results are in agreement with those previously obtained, but are significantly more 
precise. 

Once ^_L,oo(/5) has been computed, we have first performed fits of the form 

eoo(/3) = i(i-|) (68) 

using only the data with /? > /3min- The value /3mm corresponds to the lowest value for 
which the is less than the value corresponding to a 95% confidence level. Using the 
first set of extrapolated values, see Table El we obtain = 0.483(50), Pc = 0.31264(15), 
while the second set, see TableEl gives = 0.500(32) and Pc = 0.312696(88). The results 
are in full agreement with the prediction i>± = 1/2. Of course, in the second case this 
is simply a consistency check, since, by using the theoretical prediction for the scaling 
curve, we have implicitly assumed h'± = 1/2. 

In order to obtain our best estimate of Pc, we fix z/_|_ = 1/2 and use the second set of 
extrapolations (which have been performed implicitly assuming i'± = 1/2). Still fitting 
with Eq. (jUH|) . we obtain for /5min = 0.31 

/?e = 0.312694(18), (69) 

with = 6.0 and 5 degrees of freedom. The addition of an analytic correction does not 
improve substantially the result. 

The result for Pc should be compared with the existing determinations: 



Pc 



'0.3108(11) (Ref. ISni) 

0.3125(13) (Ref. ^) 

0.3156(9) (Ref. [201) 

^0.3125(10) (Ref. 



(70) 



Our result ()69|) and the estimates of Ref. j^lT^ I^ are in reasonable agreement. On the 
other hand, the estimate of Ref. is sligthly larger, although the difference (three error 
bars) is not yet very significant. 

6 Conclusions 

In this paper we have performed a thorough check of the theoretical predictions for the 
DLG. We have verified that the FSS curves F{2, z, S^) for x± ^ind agree with the 
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JSLC predictions ^ and if A is fixed to the JSLC value A = 2. We wish to stress 
that the comparison between theory and numerical data does not require any tuning of 
parameters, at variance with previous studies in which the FSS analysis required fixing 

and/or some exponents. We also analyzed the magnetization at the critical point 
finding /3±/z/^ = 1.023(43), in agreement with the JSLC value P±/i'± = 1. Our results 
for the Binder parameter g do not agree with those of Ref. where it was found g ^ 
at criticality, but confirm the results of Wang [T3] who could not find a satisfactory 
collapse for the Binder parameter. Our result = is compatible with the idea that even 
finite- volume transverse correlations are Gaussian in the scaling limit, so that g = a.t 
criticality. Finally, we discuss the FSS behavior of the magnetization. Our results are 
consistent with a Gaussian distribution for the order parameter with corrections that are 
well parametrized in terms of the Binder parameter. 

One point that is still unclear is the role of the dangerously irrelevant operator. A 
general analysis confirmed by an explicit one-loop computation indicates that this opera- 
tor does not give rise to FSS violations in any 2 < < 5. However, in two dimensions the 
same arguments predict logarithmic violations (if it is marginally relevant) or at least log- 
arithmic corrections (if it is marginally irrelevant). Numerically, we observe corrections to 
the Gaussian field-theoretical predictions but we are not able to determine their behavior 
as L±_ oo. Indeed, very large lattices are needed to distinguish logarithmic corrections 
from power-law corrections with small exponents. In any case, our data indicate that the 
operator is marginally irrelevant since g ^ as L± 0. We find '^g/i'± = 0.45(15) but 
this result should not be taken seriously. Probably, g decreases as a power of log L because 
of the marginal operator, but, in our range of values of L, the complicated logarithmic 
dependence is mimicked by a single power. Note that, if g{l3c,oo,oo) = 0, the Binder 
parameter cannot be used to compute Pc- The crossing method does not work. 

Let us now compare our results with those of Refs. [201 Ell ^^^^ presented numer- 
ical results apparently in good agreement with the RDLG scenario. If the JSLC theory 
gives the correct description of the critical behavior, the family of lattices considered in 
Ref. 1201 — fixed A = 1 — is such that L_|_ increases too fast compared to Ly. As dis- 
cussed in Ref. ^T\, in this case one expects to observe scaling with an effective geometry 
corresponding to L± = oo. Therefore, longitudinal quantities should show the correct 
behavior, i.e. with the JSLC exponents, when studied in terms of {(3 — /3c)Ly'''^" . On the 
other hand, transverse quantities should have a critical behavior with effective exponents. 
Therefore, there is no contradiction with the estimates of Ref. (201 that differ from the 
JSLC ones. It is also useful to take the opposite point of view: What should our results 
be if the DLG belongs to the same universality class of the RDLG? If this were the case, 
the correct anisotropy exponent would be A = 1. Thus, the lattices we consider are such 
that L|| increases too fast compared to L±. Again, we expect to observe scaling [211, with 
an effective geometry Ly = oo. Transverse quantities should have the correct behavior, 
i.e. transverse critical exponents should coincide with those of the theory with A = 1. 
Therefore, our results should be the same as those of Ref. which, as we have shown 
here, is not the case. Thus, our data exclude A = 1. As discussed in Ref. [HIl, it is 
possible to determine A unambiguously from the FSS behavior of and of Indeed, 
only if A is chosen correctly, both correlations length scale linearly, i.e. C,± ~ L± and 
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^11 ~ L\\. This is indeed what we have checked in Sec. I5.21 Therefore, our data strongly 
support the JSLC prediction A = 2, and do not support the claim of Ref. that the 
DLG and the RDLG belong to the same universality class. 

Finally, let us discuss |51j the results of Albano and Saracco [21]. First of all, there is 
a serious flaw in one of their scaling Ansatze, due to the fact that they do not distinguish 
between and jS]. While in their Eqs. (5), (6), and (7) the exponent z should be 
identified with z\\ (the JSLC prediction is z\\ = 4/3 as correctly appears in their Table 
1), in Eqs. (8), (9), and (10), the exponent z should be identified with z_i_ (the JSLC 
prediction is z± = 4). Therefore, their result for c± does not agree with any prediction, 
neither the JSLC one nor the RDLG one. Moreover, with the lattice sizes they consider, it 
is not clear whether they are really looking at the short-time dynamics in infinite volume 
or rather to the approach to equilibrium in a finite lattice. Indeed, since the dynamics in 
the longitudinal direction is fast, one expects to generate correlations of size Ly in a time 
of order Ly" , so that a necessary condition to avoid size effect is that t <^ " . With their 
typical lattice sizes, one would expect to avoid size effects for t < 10^. If larger values of t 
are used — they consider times up to 10^ — finite-size effects are relevant and thus, contrary 
to their claims, the aspect ratio and the value of the anisotropy exponent become again 
a crucial issue. 
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A Binder cumulant without zero mode 

In Sec. 0]we stated that the Gaussian nature of the transverse fluctuations implies that 
the Binder cumulant vanishes at criticality. We consider now the 0{N) model for N ^ oo 
above the upper critical dimension {d > 4), showing that the Binder cumulant defined in 
Eq. (fT^ also vanishes, in spite of the presence of the dangerously irrelevant operator. This 
is due to the fact that the definition (fT^ involves the correlation functions at nonvanishing 
momenta. 

To be concrete let us consider the 0{N) cr-model defined on a finite hypercubic lattice 
A of volume L"^, with Hamiltonian 

H= -NJ2<^^- ^y, (71) 

(x,y> 

where (■, ■) indicates nearest-neighbor sites. The partition function is simply 

Z= /n|d.,*(.J-l)le-^«x{l(^^^j ™ (72) 
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where f3 = 1/T, ZM corresponds to the standard case (with zero mode), and z)(d to the 
theory without zero mode. The Fourier transform of the two-point correlation function is 
given by 

G-(p)^(^>|!U^i-^x{;_,__^ ™ (73, 

where = 4 sin^(p^/2) and Xl = Xl{P,L) solves the finite-volume gap equation of 
the model 

in which A* = 2ttL~^ is the dual lattice. 

We now define two different Binder cumulants 

1 ((m ■ mf)c 



(m ■ m) 



2 



where m° = (7^^^, /i" = c'"p=p^^i„) and |pmin| = 2n/L. The definition go{P,L) is the one 
usually used in the A^-vector model, but it is unsuitable for systems without zero mode. 
In this case, the natural definition is gi{P,L) and indeed such a quantity corresponds to 
the Binder cumulant we have used in the DLG simulation, cf. Eq. (|13|). 
Now, let us compute (7o,i(/?c, L) for c? > 4 at the critical point Pc, 



In Ref. jzo] we found that, for d > 4, 

where Cd,i and X(p) are defined in Eqs. (2.19) and (A. 5) of Ref. j^E], respectively. This is 
sufficient to determine the behavior of (m ■ m)c and (/x ■ /x)c at (bulk) criticality. As far 
as the four-point function is concerned, it can be expressed in terms of 

I, ,1 1 1 J 1 ZM , , 

^ = z.^ 1(1- - w) zx. ^''^ 

After some calculation one can easily show that, at criticality and with |q| ~ 1/L, 

^ (q)U = ^(l + M + |o(L-^) zXi ('^^ 
In the large A^-limit we have 

• O-q), = A^G'l'l(q) + O(ArO), 

((^-q ■ ^^f)c = -{N + 2)A(q)[Gi'Hq)]' + 0{N'), (81) 
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so that 



9o{Pc, L) = ^^[1 + OiL'-"/', N~% ZM (82) 

g,{(3,, L) = 0{N-'L'-''). ZM, z)(d (83) 

The different scahng behavior can be traced back to the fact that, while the four-point 
function has always the same scaling behavior, irrespective of the momentum at which it 
is computed, the two-point function scales differently (for d > 4) depending on whether 
it is computed at zero or at nonzero momentum. 

It is also of interest to compute the Binder parameter in d = 4. Without reporting 
the details, we find that go{Pc, L) agrees with the result reported above, with logarithmic 
corrections. For gi{(3c,L) we obtain instead 

gML)^ ^^^'^ J" (84) 

iV^ TT^logL 

The Binder cumulant gi{(3, L) vanishes also in = 4, albeit only logarithmically. 



B Dangerously irrelevant operator and correlation 
functions 

In this Appendix we compute the zero-momentum insertion of the dangerously irrelevant 
operator A into static correlation functions of n density fields s taken at vanishing parallel 
momenta in the JSLC theory p.^ We begin by considering ip(jp,uj) = s{p\i = 0,p_L = 
p,uj), i.e. we compute 

{ip{pi,UJl) ■ ■ ■ip{pn,UJ„)A)conn- (85) 

We recall that the interaction vertex V of the theory [H] is sVys^ (s is the response field) 
and thus it vanishes whenever the parallel momenta flowing into it from the s-leg is 
zero. Let us discuss the consequences of this fact on the form of the generic diagram V 
contributing to Eq. (jHK|) . 

First, note that if the amputated diagram P'^™? associated with T> has no external 
S-legs, it contains a loop of response propagators and therefore it vanishes because of 
causality pHll34[ IH5]. Thus T> vanishes unless P^™p has at least one external s-leg. Each 
external s-leg in P'^™? is connected either to the insertion A or to an interaction vertex 
V. In the latter case, however, the contribution vanishes because of the zero parallel 
momentum flowing from the external s-leg into the vertex V. Thus, V does not vanish 
only if each of the n > 1 external s-legs of P^™p belongs to the insertion of A. Now we 
show that D^™P cannot be one-particle reducible. Indeed, let us suppose that D'^™? can be 
divided into two amputated subdiagrams X>1^^ and 1)2"^^, such that A G V^^^^, connected 
with a single line. By momentum conservation, the parallel momentum flowing into this 
line is zero and therefore such a line must be connected to an s-leg of a vertex in T>2'^^. 

^In this Appendix we use the notations of Ref. [HI ■ 
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Since we also have A G P^™^, all external legs of Dg™^ ^ire of s-type. By causality 
vanishes. 

Summing up, the amputated diagrams contributing to Eq. (j85p are one-particle irre- 
ducible and all the n > 1 external s-legs belong to the insertion of A. 

This conclusion is independent of the particular A considered. Now let us specify our 
discussion to the case of the dangerously irrelevant operator. The dangerous irrelevant 
coupling sAj^s^ mixes, under renormalization-group transformations, with other operators 
of equal or smaller scaling dimension. Those with the same dimension, that we consider 
in the following, are listed in Ref. [H]. They have the following general structure (the 
dots indicate derivatives): (a) s(---)s, (6) §{■■■)§, (c) sV||(---)s^, (d) s{---)s^. For 
n > 2, the only case we consider, the insertion of one of the operators (a) or (6) into 
the correlation function gives only one-particle reducible contributions since one of their 
s-legs has to be an external one for the amputated diagram. Thus, according to our 
previous discussion, they do not contribute to the correlation function. The same is true 
for operators (c). The s-leg, being an external one for the amputated diagram, carries a 
zero parallel momentum p\\ into the vertex. Since these vertices are proportional to p\\, 
they vanish. Thus, contributions can only come from sA^s^ G A. Defining 

B^ = -! dM'^xsA^s^ (86) 



6 . 

we obtain 

conn, amp UT i^ri-l;Bz, (87) 

where Tn^n-^Bs is the one-particle irreducible vertex function with n external s-legs, n 
external s-legs, and one (zero- momentum) insertion of the operator B^^, A ^ uB^ + . . . 
Therefore, all connected correlation functions with vanishing parallel momenta and n> 2 
are proportional to u as u — > 0. 

We wish now to compute the transverse static four-point function G^^^ (p, p, — p, — p; Ly , Lj 
at zero parallel momenta. For simplicity, we restrict the computation to r = 0, i.e. to the 
critical point. The tree- level insertion of i?3 into the four-point correlation functions with 
all fields s taken at the same value of t (we can set t = 0) and with momenta P = (0, p_l) 
is given by 



/+00 
dUR{-U,V)C\-W. 
■00 

where we have used^^ (see Ref. jH]) 



Ms 



(pi)^ 



R{ti - t2, q) = (s(t2, -q)s(tl, q))TL 

= e{ti - t2) exp{-[ql + gj](ti - t^)} (89) 
C{ti - t2, q) = (s(t2, -q)s(ti, q))TL 

= -^;;2 exp{-[qi + ql]\ti - h\}, (90) 
q_L "T yii 



^°For n — 2 (a) and (b) have nonzero tree-level contributions. 

^^We disregard dimensionless factors present in the dynamic functional of Ref. 0. 
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and 9(t) is the Heaviside step function. 

At one loop there are two diagrams whose contributions can be written in the form 

Li(ti, t2, ts; p±) = 

- E ^V^V^i^^ - ^2, q''')i?(^i - ^3, ^''')Cih - ts, q2'3), (91) 

(?ihqx) 

-^^2(^1, h, ts] p±) = 

- ^V^T^it^ - t2, q''')R{h - h, ci''')C{h - ts, q'''), (92) 

{9ii.qx) 

where Yl' runs over the whole momentum space except the zero mode (0, 0) and (0, ±p±), 
and 

q''' = (gii,q± + p±), (93) 
q''' = (-gii,-q± + p±), (94) 
q''' = (gii,q±)- (95) 

The final result is given by T = Ti + T2 where we defined 

/+00 
dtidt2dt3Li{ti,t2, h; p±)R{-h, P) Yl C{-tk, P), (96) 
fc=i 

and 71 = 4!/2, 72 = 271 are combinatorial factors. Performing the integrations over ti, it 
is easy to find 

T(px) = J2'm,p^), (97) 

with /(q, p± = 0) 7^ and /((A^gy, Aqx), p_l) ~ A^^ for large A. The sum over /(q, p±) 
is convergent by power counting (for d < 7) and thus, for \p±\ — * 0, 

E'/(q'P^)^ (98) 

(9ii.q±) (iJihqx) 

Then, if |p_|_| ~ ^/L±, it is easy to see that, for d = 5 and at S* = L\/ L\\ fixed (thus 
V = L;[+V^), T is finite and 

T ~ Ptl ~ (99) 

in agreement with Eq. ()50p. 
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Figure 1: tl = t{P; L\\, L±) for different lattices as a function of = Lj^. Filled 
(respectively empty) points refer to lattices with aspect ratio 5*2 (respectively Si) fixed. 
Here 5*2 ~ 0.200, 5*1 ~ 0.106. Errors are smaller than the size of the points. 
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Table 1: The results of our simulations with A = 2 and 5*2 ~ 0.200. Here A = 
(table continues as Table |21). 
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Table 2: The results of our simulations (continued). 
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Table 3: Integrated autocorrelation times (expressed in sweeps) for the susceptibility 
for /3 = 0.312. 
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Table 4: Fit of t(/3; -Ly, L^) with r^oiP) + a{(3)L'^ using the data with ,^2 0.200: L^nm 
is the minimum value of L used in the fit, iV is the number of degrees of freedom, R^ is 

the x^, i.e. the sum of the square residuals, and i?^ is the value of R^ corresponding to 
the 95% confidence level based on a distribution with N degrees of freedom. In the 
last column we report Ci,oo(/^) ~ l/''"oo(/3)- 
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Figure 2: Ly, L_l) for different lattices and 5*2 as a function of L^. Here a = 2, 3, 4. 
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Table 5: Results of the extrapolation with the iterative method of Ref. jlOj using the 
theoretical FSS function. Here is the number of points for each j3, B? the residual, and 
Bj^ is the value corresponding to a 95% confidence level. 
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Figure 3: FSS plot for the transverse correlation length ^± with a = 2 and 5*2 ~ 0.200. 
Here = L\\, L±), = SLy, 2L_|_), L = L^. Different symbols correspond to 
different lattices sizes: L = 16(o), 18(0), 20(A), 22(0), 24(*). The line is the theoretical 
prediction F^{2,^l/ L, S2), cf. Eq. 
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Figure 4: FSS plot for the transverse correlation length analogous to that presented in 
Fig.H Here we plot the data vs ^2l/{2L) = SLy, 2L^)/(2Lx). Symbols are defined 
as in Fig. IHl The line is the theoretical prediction 
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Figure 6: FSS plot for the magnetization m. Symbols are defined as in Fig. El The solid 
line is the theoretical prediction, cf. Eq. (jSH). 
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Figure 7: FSS plot for transverse Binder parameter g. Symbols are defined as in Fig. 01 
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Figure 8: Plot of -Ly, Li_)/L\\ versus -Ly, L±)/Lj_. Each symbol corresponds to a 
value of p: 0.28 (filled Q), 0.29 (filled □), 0.30 (filled A), 0.3025 (filled y), 0.305 (empty 
O), 0.3075 (empty □), 0.31 (empty A), 0.3105 (empty 0), 0.311 (empty cross), 0.31125 
(filled O), 0.3115 (filled □), 0.31175 (filled A), 0.312 (filled y)- For each f] points move 
leftwards as Ly and L± increase. 
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